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^ _ The effects of heavy quarks energy loss on the transverse momentum spectra of single 

, muons are studied. The energy loss of heavy quarks when traversing the medium formed 

in PbPb collisions at ^/s NN =5.5 TeV is calculated by implementing the collision geom- 
' etry and the quenching weights. The medium density is assumed to decrease at forward 

' pseudo-rapidity and to be proportional with pseudo-rapidity multiplicity dN/dn. Muons 

from W decays can be used as a medium-blind reference to quantify the effect of heavy 
i quarks energy loss on the single muon production. 

1. Introduction 

Heavy quarks are sensitive probes of the medium produced in the collision as they 
are produced in the early times and may lose energy when traversing the medium. 
At the Relativistic Heavy Ion Collider (RHIC), a significant suppression of electrons 
from heavy quarks at mid-rapidity has been ob serv ed in central AuAu collisions, 

1 1 lOI 

indicating substantial energy loss of heavy quarks^ 1 ^ 1 . At the Large Hadron Collider 
(LHC), the cm. energy of PbPb collision is 5.5 TeV per nucleon-nucleon collision, 
about 30 times larger than that at RHIC and it is expected that heavy flavors will 
be abundantly produced, offering a good opportunity to characterize the medium 
formed in the collision. Heavy quark deca ys are ex pected to dominate the lepton- 
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pair continuum up to the mass of the Z° M I 4 | o | b | / | an( j transverse momentum 
spectra of single muons up to 30 GeV/cJ§J. The shape of the spectra could be altered 
by the energy loss of heavy quarkd ^ l lO I. At the same time, this high amount of 
energy in the center of mass enables the possibility to create the W bosons in the 
initial hard scattering. As the muons from W decay are not effected by the medium 
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formed in the collision, it is interesting to see how much the energy loss of heavy 
quarks affects on the muon production by using the muons from W as a reference. 

In the present work, we study the effect of heavy quarks energy loss on the trans- 
verse momentum spectra of single muons in PbPb collisions at ,/snn = 5.5 TeV in 

[a] 

the acceptance of ALICE MUON Spectrometer . First, by introducing the nuclear 
shadowing effect and k t broadening in the HVQMNR prograrrPJ] which provides a 
fully exclusive partonic differential cross section calculation including all the rele- 
vant partonic subprocesses, we calculate the kinematic distributions of heavy quarks 
per nucleon-nucleon collision in PbPb collisions at ^/snn = 5.5 TeV. Then, to sam- 
ple the energy loss of heavy quarks, we use the collision geometry and the BDMPS 
quenching weights for massive case^l with a scaling of the transport coefficient 
as a function of the pseudo-rapidity. Finally we give the results of the transverse 
momentum spectra of single muons within the acceptance of ALICE MUON Spec- 
trometer. 



2. Production of Heavy Quarks 

In the HVQMNR program, the heavy quark hard production cross section is cal- 
culated in pQCD according to the factorization theorem: 

dcr a h(Vs) = X! / Ax ^ Ax 3 d& ^j{s,m 2 1 ^^\)f^{x l ^] ? )f b 1 {x : j^ 2 F ), (1) 

where /if is the fractorization scale, fi(xi) and fj(xj) are the parton distribution 
functions, the differential probabilities for the partons i and j to carry momentum 
fractions Xi and Xj of their respective protons a and b. The total partonic cross 
section dov, of the process ij — > QQX at parton-parton cm. energy s = XiXjS, 
where y/s is the cm. energy of the pp collision, is calculated as a perturbation series 
in cx s ((1r), where the strong coupling constant is evaluated at the renormalization 
scale fiR. 

In the present calculation, CTEQ 4M parton distributions^! are used and nu- 
clear shadowing is taken into account by using the EKS98 parameterizatiorP3l. 
The parameters of HVQMNR program are set as: for charm, m c = 1.2 GeV 
and fip = hr = 2/xq, for beauty, mi, — 4.75 GeV and fip = hr = (Iq, where 
/Jo = \J m Q + (Pt,Q q)/2- Calculations are performed up to next leading order 
(NLO) with the HVQMNR program. 



3. Energy Loss of Heavy Quarks 

For modelling the energy loss of heavy quarks, we used the quenchi ng w eights in the 
multiple soft scattering approximation, which were derived in RefM^ 1 in the frame- 
work of the BDMPS formalisrrPS] The energy loss AE is sampled by i mp lementing 
the collision geometry and the quenching weights as explained in Ref. . 
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However, the transport coefficient q mentioned in the RefP 

is only applied for 

the central rapidity region. The decreased medium density in the forward rapidity 
region should be taken into account. We assume the scaling 
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where dN^/drj is the pseudo-rapidity distribution of charge particles predicted 
17 

in Ref. . The transverse momentum of the quark in the medium is reduced to 
p T = pt — AE, while the rapidity is kept unchanged. 

Heavy quarks that lose most of their energy due to medium effects are redis- 
tributed according to a thermal distribution^, 

dA^thermal 



TOT . 



where we choose T 



dTOT 

300 MevEJ 



oc m T exp(--—) 



(3) 



4. Heavy quarks contribution to single muons 

The ha dron ization of heavy quarks is assumed to follow the Peterson fragmentation 
function^, 

Dh ^ {z)cC z[l-l/z-\ Q /(l-zW (4) 

where z = Ph/pQ is the fraction of heavy quark momentum carried by the final 
state hadron. eg is a parameter f or t he fragmentation, e c = 0.021 for charm and 
£b = 0.001 for beauty are used here'I^. 

Heavy mesons produced after the fragmentation decay semileptonically. In the 
spectator model, the heavy quark in a meson is considered to be independent of 
the light quark and allowed to decay as a free particle with the simple V-A weak 
curren Thus, we considered heavy quark three-body-decay(c — > /i + v + s and 
b — > fi+v + c) into muons instead of heavy meson decays into muons. The decays are 
calculated in the rest frame of the heavy quark and then boosted back accor ding to 
the heavy quark momentum to obtain the kinematic distributions of muon; 
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5. Results 

Fig. Q] shows the transverse momentum distribution of single muons from heavy 
quarks in pp collisions at ^/snn = 5.5 TeV obtained from PYTHIA simula- 
tions (PYTHI A has been tuned to reproduce the kinematic distribution of heavy 
quarks produced by the HVQMNR program with parameter setting described 
in Section 2) and our Fast Simulation composed of the HVQMNR program, Pe- 
terson fragmentation and the spectator model(mentioned in the Section |4|). k t 
broadening and nuclear shadowing are included. In Fig. [T] the left plot is with- 
out any cuts and the right one is within the acceptance of the ALICE Muon 
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Fig. 1. The transverse momentum spectra of single muons from heavy quarks in pp collisions at 
y/sNN = 5.5 TeV obtained from the PYTHIA simulation and the Fast Simulation composed of the 
HVQMNR program, Peterson fragmentation and the spectator model, kt broadening and nuclear 
shadowing effects are included. The left plot is with no cut, the right one is within the acceptance 
of ALICE Muon Spectrometer. Muons from the decay B — » D — > /i are not taken into account. 
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; q(0) = 

q{0) = 25 GeV J /fm, with dN/dTi scaling 

q{0) = 100 GeV 2 /fm, with dN/dii scaling 

q{0) = 25 GeV J /fm, without dN/dri scaling 

; 4 q(0) = 100 GeV 2 /fm, without dN/dTi scaling 
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Fig. 2. Nuclear Modification factor for heavy quarks in the central(0-5%) PbPb collisions at 
\/ s NN = 5-5 TeV. Lines for different q values with/without dN/dn scaling are given. The left plot 
is for beauty and the right one is for charm. 



Spectrometer(2.5 < t/ 1 < 4, > 4 GeV/c, p£ > 1 GeV/c). Muons fr om the 
decay B — > D — > fi are not included. The Fast Simulation gives transverse momen- 
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Fig. 3. Transverse momentum spectra of single muons in central(0-5%) PbPb collisions at ^/snn = 
5.5 TeV within the acceptance of ALICE Muon Spectrometer. Transport coefficients are all applied 
with dN/dn scaling. The pt spectra of all muons and the muons from W are obtained from RcfA 
The pt spectra of the muons from heavy quarks are obtained by the Fast Simulation. 



turn distributions of single muons from heavy quarks close to those obtained from 

pythiaP. 

To see how much the transport coefficient dN/drj scaling (formula (2)) affects 
the energy loss of heavy quarks, nuclear modification factors for different values 
of q with/ without dN/d-q scaling are shown in Fig. [2] The left plot and right plot 
are for beauty and charm, respectively. From Fig. [2l we can get that the decreased 
energy density in the forward rapidity region makes very small effects on the energy 
loss of heavy quarks. 

Fig. [3] presents the transverse momentum spectra of single muons in central(0- 
5%) PbPb collisions at ^/snn = 5.5 TeV within the acceptance of ALICE Muon 
Spectrometer(2.5 < rf < 4, p^ > 4 GeV/c, > 1 GeV/c). The p t spectra of 
all muons and the muons from W are obtained from Refill The pt spectra of the 
muons from heavy quarks are obtain ed by the Fast Simulation. The values of q(0) 
are set to be 25, 100 GeV 2 /fi£2UM. As the muons from W bosons are not effected 
by the medium formed in the collision, they can be used as a reference to quantify 
the effect of heavy quark energy loss on the single muon production. From Fig. [3l 
we can see that when energy loss is applied to heavy quarks, the crossing points of 
the pt spectra of single muons from heavy quarks and from W are shifted to the 
lower ^(when q(0) — 100 GeV 2 /fm, for beauty, the crossing point is shifted from 
23 to about 17 GeV/c, for charm, it is shifted from 18 to about 12 GeV/c). 
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6. Conclusions 

We obtain similar pt spectra of single muons from a Fast Simulation(composed 
of the HVQMNR program, Peterson fragmentation and the spectator model) and 
from PYTHIA generator. Accounting for the decrease of the medium density in 
the forward rapidity region, we assume a dN/dr] scaling to q and we find that the 
decreased energy density in the forward rapidity region makes small effects on the 
energy loss of heavy quarks. The energy loss of heavy quarks causes the crossing 
points in p t spectra of single muons from heavy quarks and from W bosons shift to 
lower p t . 

Acknowledgements 

The authors would like to thank F. Antinori and the ALICE group members in 
Padova/Legnaro, Italy. And G. Martinez-Garcfa is gratefully acknowledged for 
his helpful discussion and enthusiastic help for the work. This work is supported 
partly by the Grant of Museo Storico della Fisica e Centro Studi e Ricerche 'En- 
rico Fermi', the National Natural Science Foundation of China under Grant No. 
10575044 and Key Grant No. 10635020, the Key Project of the Chinese Ministry 
of Education No. 306022, and partly by the EU Integrated Infrastructure Initiative 
Hadron Physics Project under contract number RII3-CT-2004-506078. 



References 

1. S. S. Alder et al, Phys. Rev. Lett 96 (2006) 032301 

2. X. Dong et al., Nucl. Phys. A774 (2006) 343. 

3. F. Carminati et al., J. Phys. G: Nucl. Part. Phys. 30 (2004) 1517. 

4. B. A. Kniehl, K. Originos and I. Sarcevic Nucl. Phys. A700 (2002) 523. 

5. I. P. Lokhtin and A. M. Snigirev, Eur. Phys. J. C21 (2001) 155. 

6. D. C. Zhou and D. Jouan, IPNO-DRE. 96-24. 

7. D. C. Zhou, et al., Chin. Phys. Lett. Vol.16, No.7 (1999) 496. 

8. Z. Conesa del Valle et al., ALICE-INT-2006-021. 

9. Z. W. Lin, R. Vogt and X. N. Wang, Phys. Rev. C57 (1998) 899. 

10. H.-T. Ding, D. C. Zhou and Z. B. Yin, Chin. Phys. Lett. Vol.23,No.9 (2006) 2399. 

11. M. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B373 (1992) 295. 

12. N. Armesto et al., Phys. Rev. D71 (2005) 054027. 

13. H. L. Lai et al., Phys. Rev. D55 (1997) 1280. 

14. K. Eskola, V. Kolhinen and C.A. Salgado, Eur. Phys. J. C9 (1999) 61. 

15. R. Baier, et al., Nucl. Phys. B483 (1997) 291. 

16. A. Dainese, C. Loizides and G. Paic, Eur. Phys. J. C38 (2005) 461. 

17. D. Kharzeev, E. Levin and M. Nardi, Nucl.Phys. A747 (2005) 609. 

18. C. Peterson et al., Phys. Rev. D27 (1983) 105. 



19. J. Baines et al., hep-ph/0601104 



20. G. Altarelli et al., Nucl. Phys. B2 08 (1982) 365. 

21. T. Sjostrand et al., |hep^ph/0108264| 



22. N. Carrer and A. Dainese, hep-ph/0311225 



